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A standard two-path interferometer fed into a linear Al- 
port analyzer with coincidence detection of its output ports 
is analyzed. The AT-port is assumed to be implemented as a 
discrete Fourier transformation T, i.e., to be balanced. For 
unbound bosons it allows us to detect JV-particle interference 
patterns with an JV-fold reduction of the observed de Broglie 
wavelength, perfect visibility and minimal noise. Because the 
scheme involves heavy filtering a lot of the signal is lost, yet, 
it is surprisingly robust against common experimental imper- 
fections, and can be implemented with current technology. 
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A quantum state containing N particles can yield an 
AT- fold reduction of the observed de Broglie wavelength in 
an interference experiment when compared with a single 
particle pattern [1]. This can lead to an A^-fold increase 
of interferometric sensitivity, in the best case, reaching 
the quantum mechanical 'Heisenberg- limit' of particle- 
number-based enhancement in resolution [2] . For bound 
particles this is straightforward [3], as long as they are 
not split up by the interferometer, namely, their binding 
energy exceeds the interaction energies encountered when 
passing the interferometer's beam splitters and merg- 
ers [1]. But, it is difficult to reach the Heisenberg-limit, 
or, more generally speaking, to detect A^-particle pat- 
terns in imaging or interferometry using unbound parti- 
cles [4-6] such as free photons. 

Experimental results have so far only demonstrated 
a halving of the observed wavelength using two-photon 
states [7-9] or reported signatures of four-particle en- 
tanglement for photons [10,11] and ions [12]. Yet, 
weakly or unbound bosons such as photons or cold 
atoms are currently the most important particles for 
interferometry and imaging. Interest in several recent 
schemes to employ their multi-particle features in inter- 
ferometers [1,4,5,7-9,13], high resolution imaging [14], 
and quantum lithography [15] has therefore been con- 
siderable [16]. 

Observing multi-particle quantum effects for unbound 
particles is hampered by several problems: special 
quantum states are needed that are hard to synthe- 
size [1,2,5,14,15], the processing involves large non-linear 
particle-particle interactions [1] or joint detection of all 
particles with single particle resolution [5,7-9,13,15]. 

The linear scheme presented in this paper shows that 
equipping an interferometer's output with a balanced Al- 
port analyzer allows us to circumvent several of the above 
problems and detect multi-photon interference patterns 



with current technology. Note, that other linear schemes 
have recently been devised to circumvent similar prob- 
lems in quantum information theory [17,18] and quantum 
state preparation [19,20]. Using the scheme presented 
here it should be straightforward to observe, for the first 
time, a more than two-fold reduction of the effective de 
Broglie wavelength of unbound photons [1,7-9,13]; the 
scheme is sketched in FIG. 1. 
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FIG. 1. Sketch of the setup: dotted lines depict N — 2 
empty (vacuum) modes fed into the Fourier mixer Tn to- 
gether with the two active interferometer modes a\ and a\ 
resulting from mixing the modes a' and (3' at a balanced 
beam splitter. Mode a\ suffers the interferometric phase shift 
<j>. The corresponding output modes of the Fourier mixer Tn 
are A\ to ^4^; they are detected by N single-photon detectors 
D which are read out in coincidence. 

In order to model our system let us follow the evolu- 
tion of a state of bosonic particles starting out in chan- 
nel a' of FIG. 1 and progressing from left to right. For 
the transformations describing the mixing of the inter- 
ferometer's input modes by a balanced beam-splitter in 

conjunction with the phase shift A r L " 

the following operator equations 



for mode a\ we choose 



a\ = (&* + /? f )/V2, 4 = e~ # (a t -/3 t )/V2- (1) 

Subsequent to the interferometer the balanced Alport 
Tm acts as a discrete Fourier-transformer [21,22] and 
mixes the interferometric modes a\ and a\ with the N — 2 
vacuum modes aj, . . . , ajy . Using a vector notation for 
these input modes a\ and the output modes At, i.e. 

A' = J-n &< j we specify the corresponding (Fourier-) 
transformation matrix elements as 

1 r .27T 



{TN)lk = 7N ^pt^'- 1 )^- 1 )] 



(2) 
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To make a connection with the textbook case of classical 
or single photon interference [23] let us look at the case 
N = 2, the conventional Mach-Zehnder interferometer. 
For the input state we assume a single photon is enter- 
ing through channel a\ namely |V>)in = &t|0). After 
the transformations (1) and (2) this state is converted 
into \ip) ln = (1 + e"^)/2 (A\ - i|)|0) yielding the cus- 
tomary classical interference pattern (Ii) = (A\Ai) = 
(1 + cos</>)/2 and the antifringe pattern in the second 
channel (I2) = (1 — cos</>)/2 . Of course, the coincidence 
pattern (Ii • 7 2 ) =0 vanishes, since only a single pho- 
ton is present. The simplest non-trivial case is that of 
the Mach-Zehnder interferometer fed with a two-photon 

state \ip(2))in = ^=-|0): it yields the classical [23] first 
order interference patterns (I1/2} = 1 ± cos</> and the co- 
incidence signal (Ii -12) = (1 — cos 2^) /4, which shows 
the desired halving of the effective de Broglie wavelength 
and full contrast [23] observed in experiments [7-9]. 

Generalizing the above discussion we now consider the 
N- channel coincidence count operator In 
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where we have used the fact the the output modes A 
commute. Expressing flfcLi ^fe m terms of the input op- 
erators a\ and a\ yields 
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(4) 
(5) 



where we have used (2) and assumed that a\ = a\ = ... = 
a N = in accordance with our scheme's stipulation to 
leave these modes empty (tracing out the vacuum state) . 

Note, that the terms "a" represent generalized N- 
particle intensity or TV-particle dosage terms [15]. Cor- 
respondingly, the terms in square brackets represent 
./V-particle cross-mode terms which generalize the well 
known single particle (or classical) interference terms [23] 
of the form a\a 2 + aia\. 

Eq. (5) contains a good and a bad message, the bad 
aspect is the emergence of the signal suppression factor 
N~ N which is due to the fact that we require all output 
detectors to fire - a rare event: (In) = if one of the 
output channel detectors fails to fire. Such losses are 
common whenever one tries to substitute non-linear by 
linear elements [17-20]. 

The positive aspect is the result that the form of In 
in Eq. (5) suggests that for those subensembles of events 
that trigger firing of all N detectors a 'perfect' N^ 1 or- 
der interference pattern measurement can be performed. 
Only linear elements and an array of conventional de- 
tectors is needed for its implementation, it will turn out 
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that the scheme is rather robust with regards to detec- 
tor imperfections. Current technology suffices, Eq. (5) 
therefore is important, it is the main result of this paper. 

In order to prove that In does indeed describe a 'per- 
fect' iV^h order interference pattern measurement let us 
check its signal and noise properties. The combination 
of operators suggests that In probes for the presence of 
the number-entangled state 
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This is the superposition-state needed when trying to 
reach the Heisenbcrg limit in interferometry [2] using sin- 
gle photon count techniques. Note, that to synthesize 
this kind of state one typically needs correlated input in 
both modes a} and /r, see e.g. [20], this case is obviously 
not covered by our sketch of the setup in FIG. 1 where 
the source is only seen to feed into mode 6^ . 

Its signal is (I N ) = N\/N N x (1 - (-l) N cos(N<f>)) and 
does indeed display the term cos(iV</>) due to an iV-fold 
decreased effective de Broglie wavelength [1,7-9,13] and 
an interference pattern with perfect visibility. 

The second moment is (1%) = N\ 2 /N 2N sin 2 (N<j)). 

To study the scheme's noise features, we will estimate 
the noise induced phase spread A(j) from the quotient of 
the signal's fluctuations and the signal's phase gradient 
(see a good textbook such as [24]) 



A<j> = M N /\dX N /dcj)\ 



(7) 



with A2n = 



Xn) ~ @n) 



For the number- entangled 



state (6) the phase spread therefore is 
Ac/) = V-W • 



(8) 



The subensemble of registered events (all detectors fire) 
reaches the Heisenberg-limit [24] which proves that our 
scheme has the least possible noise [2]. 

According to Ou's analysis [2] this limit cannot be 
reached using arbitrary states, i.e. in our case, states 
other than (6). 

Let us therefore study the signal and noise properties of 
some representative states that are fed through one chan- 
nel, e.g. a\ only. To start out, let us consider iV-photon 

Fock states \ip(N))in = ^7=; |0), which can in principle 
be synthesized using projection measurements [19]; the 
signal (In) has the form 
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&n) = (U^ = 2 (n-i) N n (l-(-l) N coa(N<f>)) . (9) 



We see that an interference pattern (9) with TV-fold re- 
duction and perfect contrast is observable, yet we will 
find that the imperfect matching to the ideal state (6) 
leads to large noise. 
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According to Stirling's formula, this signal scales like 
with the corresponding numeri- 
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cal values 1/4, 1/18, 3/256, 3/1250, 5/10368 for N = 
2, 3,..., 6. Evidently, the requirement for simultaneous 
firing of all detectors reduces the signal strength consid- 
erably. 

The joint-detection operator In extracts suitable com- 
ponents to show iV-photon interference patterns out of 
every state with sufficiently large photon numbers. In 
order to prove this for general input into channel cr let 
us now consider an excess E of particles N sta te = N + E 
over the number of ports N . Trivially, a negative ex- 
cess, i.e. a number of particles fewer than expected de- 
tector clicks will give a zero signal. For a Fock-state 

= \ ' 10) with a positive excess the cor- 

responding intensity expression shows the same pattern 
as Eq. (9) above - with favorably increased magnitude 



(N + E)\ 1 - (-1) N cos(N(f>) 
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In order to appreciate the very considerable intensity gain 
attributable to such a photon-number excess, note, that 
at the rather modest cost of doubling the number of pho- 
tons, one can roughly compensate the losses due to an in- 
creased number of ports. For example, the losses due to 
moving from N = 2 photons and channels to 10 photons 
and channels is offset by moving to 10 channels and 20 
photons instead. Moreover, large photon number excess 
also reduces the noise, see below. 

In order to generalize our results to states other than 
Fock-states let us maintain that the /3-channel is empty 
and let us trace it out: (3 = ft = 0. In this case 
the coincidence count operator In commutes with the 
input photon number operator era implying that any 
Fock-state component exceeding the port number N con- 
tributes separately to the overall interference pattern and 
its noise features. Hence, for a general a-channel input 
state \ip)i n = SjLq Cj ^Jf |0)> * ne coincidence pattern 
has the form 



{I N ) = \cn+e\ 



E=0 
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A corresponding generalization applies to mixed input 
states; this shows that our scheme allows for a robust 
formation of an A-fold reduced interference pattern with 
perfect visibility in the general case of states fed through 
channel a) only. 

Let us now discuss the noise properties of the N- 

photon state \ip(N))in — ^=p|0). The second moment is 

(iff) = M 2 /(N 2N 2 N - 2 ) x(l - (-1) N cos(N<t>)) yielding 
the minimum phase estimate (at positions = for even 
N and 4> = tt/N for odd N) 



Acf) = 



N 



(12) 



A Fock state fed into the a^-channel has little overlap 
with the number-entangled state (6). This leads to extra 
noise increasing the phase spread by the factor ■\/2 N ~ 1 . 

We have seen that the use of Fock states with a photon 
number excess N + E over the port number N leads to 
a considerable signal increase (10). Also the noise per- 
formance in this case is better. The corresponding ex- 
pression for the noise is a bit involved (it was determined 
using the 'Maple' symbolic manipulation program). For 
large photon number excess E we find that the noise falls 
to the shot-noise level, but not below; see FIG. 2. 
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FIG. 2. Logarithmic plot of minimum phase spread Acf) 
(logarithm 'Id' to basis 10) for setups with N ports. The 
poor performance (12) of a Fock state with N photons fed 
into the aL cnanne ] (top sheet) is compared to shot noise 
Acf) — 1/%/iV + E (bottom sheet). Shot noise can be reached 
employing coherent states \yN + E) a . For a Fock-state 
\ip(N+B))in with sufficiently large photon excess E over the 
number of ports N (middle sheet) the noise is reduced, reach- 
ing (but not falling below) shot noise. 

Note, that we restricted most of the above analysis to 
the case of states fed through the input channel aJ only, 
this is because there are few transparent analytical ex- 
pressions in the general case. However, in order to go 
below the shot noise level quantum correlated input, us- 
ing channel ft as well, is needed (remember Eq. (8)); 
this can either reduce or enhance the interferometric sig- 
nal [23]. Our discussion shows that currently available 
quantum states can be used and that moving towards the 
two-mode entangled number state [20] allows us not only 
to observe an A-fold reduced de Broglie wavelength but 
also noise below the shot-noise level (for the subensem- 
bles selected by the coincidence clicks). 

Our scheme is not only versatile with regards to the 
states that can be used, it also is quite insensitive to 
detector imperfections (other than dark counts) which 
often lead to serious degradation for quantum state re- 
construction schemes [25]. Let us first consider detec- 
tor losses, described by the admixture of vacuum modes 



V, in* every detector mode A 



t 

1. We find that the multi-channel coinci- 



T* A 1 
T 3 A 3 
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dence signal operator for lossy detectors In, lossy acquires 
the form 
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N 

±N,lossy = l[(r* A] + pltyfaAj + Pj V S ) (13) 
3 = 1 

N N N 

= n n A J A i = t n ^ = . (w) 
fe=i j=i j=i 

Again, the fact that all modes Vj are vacuum modes 
and can be traced out by setting Vj = has been used. 
The only net-effect of detector losses is the rescaling of 
the signal strength by the transmission loss factor T = 
rifeLi \ T k\ 2 ■ Since this factor cancels when forming the 
appropriate quotients of Eq. (7) it does not even affect 
the phase-resolution properties directly. 

Most detectors are unable to distinguish one from two 
and more photons, see [26] though. This inability of- 
ten poses problems [5,24,25] which our present scheme 
is also quite insensitive to. This can be seen from 
the following argument: if only one photon exits per 
channel Aj there obviously is no difference in detec- 
tion by number-sensitive detectors and those that only 
discriminate between presence and absence of photons. 
For the later case let us, in analogy to expression (3) 
define the N-channel photon-presence operator pv = 
Il^Li (lj — wner e the projectors 'unity minus 

vacuum' (lj — |0)jj(0|) project out the vacuum modes 

in the exit channels Aj . This is equivalent to saying that 
we determine the signal of the ./V-channel coincidence 
count operator In divided by the respective photon 
number in the output pv = Pv^JvdljLi A ] A j)~ 1 pN- 
Consequently (Pv(</>)) < @n {<!>)) and (P N (<f>)) > iff 
(2"jv(</))) > 0. In other words (Pv((/>)) and (2jv(^)) have 
the same zeros and similar global behavior: the same pe- 
riod in interference patterns can be seen. It is not difficult 
to show that the diminished signal (pv (</>)) encounters 
increased noise though. 

Note the superficial similarity of our setup to imaging 
in the Fraunhofer-limit, see FIG. 1. The 'object' (beam 
splitter plus phase shifter) is being mapped by a 'lens' 
(TV-port) into the 'far-field' (detector array Aj). It re- 
mains to be seen whether the discussion presented here 
can be extended to quantum enhanced imaging [14]. 

It should also be interesting to investigate our setup 
for correlated input states feeding several of the input 
modes. 

To conclude: we have studied an interferometric setup 
using an iV-port readout containing only linear elements. 
The N output detectors have to fire in coincidence, this 
way a lot of signal and noise are filtered away. The se- 
lected subensemble can reach the Hciscnberg limit which 
proves that the scheme is noise-optimized. The scheme 
is versatile with respect to the type of quantum states 
that can be employed, allows us to see an A-fold reduc- 
tion of the measured de Broglie wavelength with perfect 
visibility, is surprisingly insensitive to common detector 
imperfections, and readily implemcntable using present 



technology. 
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